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A DIRECT METHOD FOR THE MEASUREMENT 
OF CONTROL ROD WORTH IN A NUCLEAR REACTOR 


By 


William Arthur Emslie 


It is often necessary to perform control rod calibration 
in a nuclear reactor. Such calibration requires the measurement 
of the reactivity produced when discrete sections of a control rod 
are withdrawn from a critical reactor pile. From these data a 
total reactivity worth curve can be complied for the length of the 
rod. The subject of this thesis is the development of a means by 
which reactivity worth can be measured. 

A solution to this measurement problem has been developed 
using the current output of a compensated ionization chamber 
(adjacent to the reactor) as a signal source. When a section of 
control rod is withdrawn from a critical reactor, the current from 
the ion chamber is 

t/9 


: да р 
i (t) те 


where t is time, өв is the reactor stable period and Io is 


the initial ion chamber current. 


The stable period is proportional to the amount of reactivity 


produced since 
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9 = So 
Р ka 


Оо | D 
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9% is the mean neutron lifetime, к is the excess multiplication 


factor and p is the reactivity. This approximate relationship 


between 9 and p makes it possible to obtain the reactivity 


from the output current of the ion chamber since 


t t 
ДА р 18, 


i (t) = Ig? = Ig? 


Actual measurement of р is accomplished with a circuit 
composed of a logarithmic amplifier followed by a differentiator 


and filter. The transfer function of this circuit is 


(-R4C,s) 1 


V (в) = V. (o) me e oe 
out 1 (R4C st 1) (RC s t1) (В С 531) 


where V (в) is the frequency domain expression for У 160) (the 


output voltage of the logarithmic amplifier) and is equal to 


_ -А log, 5 (Ig/I,) А log, ge 

V4 (s) b снн ee O 
S 05 
Р 


In the time domain, the logarithm of the ion chamber current 


produces a linear ramp voltage v1 CO where 


i (t) Іе 
Т 778198071 


Ж 
Г k - 


v(t) = -A 19810 = -A Іов ут - à, log, ge 


қ 


This ramp voltage is then differentiated and filtered to 
obtain a d-c voltage proportional to the inverse of p. The re- 


sulting voltage (V5) is 
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I 
d d 0 At 
ee Сы aay LEE ыш. А 
Var GOD E CK А оо en 
k p 
_ КА pKA 
o log, ye = 8 log, ge 


V, may be calibrated on a d-c voltmeter to indicate the 
corresponding reactivity. 

À circuit to synthesize the transfer function was designed 
and constructed. Results of tests made on circuit performance show 


the total accuracy in the measurement of p to be within + 0.5 


percent of the meter scale. 
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СНАРТЕК І 


INTRODUCTION 


The necessity to perform calibration of the control rods 
in a nuclear reactor requires a means by which to measure the 
reactivity produced by withdrawal of discrete sections of a con- 
trol rod from a critical reactor pile. 

Use of an ionization chamber adjacent to the reactor pro- 
vides an exponential current (i (t)) from which the reactivity 
may be derived. When a section of control rod is withdrawn from 
a critical reactor i (t) becomes 

t/Q 


E Р (1) 


where Ig is the initial ion chamber current, t is time and 9 
P 


e 
is the reactor stable period. 
The reactivity may be obtained from equation (1.1) since 


Ө 

О 
er 1% 
95 75 (1.2) 


where 99 is the mean neutron lifetime. 


1 ) 
The reactor stable period used in this expression is the time 


required for the ion chamber current (i. Ct) to increase by a 
factor of e. 


This approximation is explained and justified in Section 2.1. 





The relationship between p and i (t) expressed in 
Equations (1.1) and (1.2) makes it possible to obtain р from 
i (t). This may be accomplished by means of a circuit which first 
takes the logarithm of i (t) and then differentiates and filters 
the result to yield a d-c voltage proportional to p. 

Àn example of a control rod calibration curve is illustrated 
in Figure 1.1. The curve is compiled by withdrawing the rod a 
short distance and measuring the corresponding reactivity produced. 
This procedure is repeated until the entire length of the rod has 
been EU 

Chapter II of this study details the theory underlying 
reactivity and the ionization chamber. The effects of prompt and 
delayed neutrons are also discussed. 

In Chapter III the reactivity measurement problem is 
formulated with a set of possible solutions. А particular solution 
method from these alternatives is then presented. Analysis of the 
design, mathematical properties, and calibration of the circuit follows. 
Further design analysis is presented in the Appendices. 


3 


The units of reactivity as expressed in Figure 1.1 are in dollars 
and cents where 1 dollar - 100 cents and is defined as 


1 dollar = XL (1-3) 
В; Bi 


where k x is the reactor excess multiplication factor and В, 15 


a normalizing constant equal to the total fraction of delayed neutrons 
and varies with different reactors. For the Michigan State University 
TRIGA Mark 1 reactor 8.+ 0.0073 and therefore one dollar's worth 

of reactivity is that amount required to produce an excess 
multiplication factor of 0.0073. 





p, Reactivity (dollars) 


4.00 


3.00 


“00 


1.00 


ddl. p ue o ls Uc. occ ee Er SE 
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Rod Position (arbitrary units) 


Figure 1.1 Calibration curve for the shim rod in the Michigan 
State University TRIGA Mark I nuclear reactor 





Chapter IV describes a prototype instrument which was built 
from the design in the previous chapter. Test results and circuit 
specifications are also given. 


Conclusions are made in Chapter V. 





СНАРТЕК ІІ 


REACTOR TRANSIENT BEHAVIOR 


2.1 Reactivity, Reactor Period and Theoretical Equations 


The excess reactivity (hereafter referred to as reactivity) 
in a nuclear reactor is defined as the ratio of excess to effective 
multiplication factors and is given by 


x 
о = - (2.1) 
eff 





In the above formula the effective multiplication factor, 
k cf is equal to unity for a critical reactor. The excess multi- 


plication factor is 
k =k - 1 (222) 


and is equal to zero for a critical reactor. 


Substitution of Equation (2.2) into Equation (2.1) yields 


k ef к 
oa 


(275) 


Since the values of k used in this analysis are to be no 
e 


ЕЕ 
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The content of this section, except for that material noted 
otherwise, is taken from M.M. El-Wakil, Nuclear Power Engineering 
(New York: McGraw-Hill Book Co., 1962), pp. 128-130. 





greater than 1.037, Equation (2.3) may be approximated as 
р =k PS (2.4) 


Using the above approximation for reactivity it is possible 
to relate р to the instantaneous neutron flux within the reactor. 
The equation describing the flux in a supercritical reactor 


is 
ф = ф е р (2.5) 


The reactor stable period (65) may be written as 


Ө 
— (2.6) 
р ех 
and substitution of Equation (2.4) for k = relates р to > 
е 
where 
Ө 
ao ee (1.2) 
P P 


This approximate expression for Е may be substituted 
into Equation (2.5) to form a relationship between reactor flux 


and reactivity such that 


t 
р ДА 


© = 556 (227) 


Reference to the rod calibration curve in Figure 1.1 shows the 
maximum reactivity to be approximately four dollars. This value 
corresponds to an effective multiplication factor of 1.0292 as 


defined in the Nuclear Reactor Operations and Training Manuel, 
Reactor Theory (Michigan State University, 1969), pp. 4-5. 





2.2 Relationship Between Reactivity and the Reactor Stable Period 


As Equation (1.2) illustrates, the approximate relation- 

ship between the reactor stable period and the reactivity is weighted 
by the mean neutron lifetime. This equation is characterized by a 
nonlinear property of 99 which is not constant but rather varies 
slightly with p. Examination of the exact expression for reactivity 
reveals the source of this nonlinearity. р may be written exactly as 

EU + 2 M , 

р да. Dit (2.8) 

where В. is the fraction of delayed fission neutrons in the ith 
group and he is the respective decay constant. 


For large 95 Equation (2.8) may be approximated чай 


ік 
ла?) 
о 

ED) 


(20) 


ср 


р = ee 
pot 


This approximation, however, is only valid for 2 > 100 and as 
EE decreases in value both the first term on the right in Equa- 
tion (2.8) and the 1 in the denominator of the summation no longer 
remain negligible. The net effect is a nonlinear relationship 
between 9, and p as illustrated in Figure 2.1. 


2.3 Effects of Prompt and Delayed Neutrons on Reactor Transient 
Behavior 


The fission process in the Michigan State University TRIGA 


Mark I nuclear reactor is described by the following nuclear 


Samuel Glasstone, Nuclear Reactor Engineering (New York: 
D. Van Nostrand Co. Inc., 1963), pp. 244-245. 


Ibid., pp. 244-245. 





p, Reactivity (cents) 
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E Period (seconds) 
Figure 2.1 Relationship between reactivity and reactor stable 


period for the Michigan State University TRIGA Mark I 
nuclear reactor (Courtesy of Gulf General Atomic) 


1000 


o 





equations: 


235 1 236 
— 2.1 
929 + o? 920 + energy (2.10) 


от 


je? + оп. - 2 fips + energy + y + ж (2-11) 


Equation (2.10) occurs about twenty percent of the time. The 
fission products (fips) in Equation (2.11) have a variety of atomic 
weights but will generally weigh about half as much as the ae 
Heat is generated by the energy given off in both reactions and the 
energy is about 200 Mev per fission. On the average between two and 
three (2 < с s 3) neutrons are given off per reaction in Equation 
(2.11). This equation is critical to reactor operation since the 
neturons produced must sustain the fission chain reaction. 

Neutrons emitted in Equation (2.11) are of two types: 
Prompt neutrons and delayed neutrons. Prompt neutrons make up 
about 99.27 percent of all neutrons emitted in the TRIGA reactor. 
The total mean lifetime of these prompt neutrons is less than one 
As e condi” 

Delayed neutrons are the result of the decay of certain 
fission products in Equation (2.11). An example of a delayed 


neutron emission is given by the nuclear equation 


3 
nu Der I р = ate” (2.12) 
22 сес Е instantaneous е 


Training Manuel, Reactor Theory, p. 1. 


Ibid., p. &. 
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where ыты is one of the fission products. The 22 second delay 


in the neutron emission gives it a half life of appearance of 22 
endow 

Even though the fraction of delayed neutrons is small, their 
long apparent lifetimes make the average lifetime of all neutrons 
emitted much longer than the lifetime of prompt neutrons alone. 

This effect is vital since a reactor period caused by prompt neutrons 
alone would make control of the fission process impossible. 

The effect of prompt and delayed neutrons on the change in 
neutron flux can be described as follows: In a critical reactor 
withdrawal of a control rod some discrete distance will produce a 
transient period due to prompt neutrons and a stable period due to 
the combined effect of prompt and delayed neutrons. The effect of 
the transient period will result in a "bulge" in the neutron flux. 


This effect is shown graphically in Figure DoD 


p positive 


іш 


p negative 


time (seconds) 


Figure 2.2 Combined effect of prompt and delayed neutrons in 
a supercritical reactor 


= El-Wakil, Nuclear Power Engineering, p. 131. 
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Since p is proportional to the slope of the curves in 
Figure 2.1 it is necessary to let the transient period die out 


prior to measurement of р. 


2.4 The Compensated Ion Chamber 


It has been stated in Equation (2.7) that the reactivity 
produced when a discrete portion of control rod is withdrawn from 
a critical reactor pile is exponentially proportional to the change 
in neutron flux. In order to measure p it is therefore necessary 
to monitor this flux. This is accomplished in the TRIGA Mark I 
reactor by means of the output current of a compensated ion chamber 
located near the reactor core. 

The output of the ion chamber is a very low current which 
is directly proportional to the neutron flux in Equation (2.5) 
such that 


an 
1 t = І Ф 
i (t) = Le aaa 
where 9 may be substituted by the approximate expression 
Ө 


ө = 2 (1.2) 
as previously stated. 

Figure 2.3 is a simplified diagram of the compensated ion 
chamber. The device is actually composed of two chambers where 
Ro is the equivalent output resistance. One chamber is lined 
with boron enriched with p.» and is biased at + 580 volts. From 


this sub-chamber a current proportional to neutron flux and gamma 


flux is obtained. The other sub-chamber has a variable bias from 
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O to -37 volts and is sensitive to gamma flux only. If the gamma 
sensitive current is subtracted from the gamma, neutron current 
of the boron coated sub-chamber, a current proportional to neutron 


flux alone is re 





Boron 
Coating 


Equivalent Output 
Resistance 





Figure 2.3 Simplified diagram of a compensated ion chamber 


At low neutron flux levels (6 « 2.5 X jouer ом the 
current contribution from the gamma sensitive chamber is nearly 
equal to that of the gamma, neutron sensitive chamber. It is 
therefore important to set the compensating chamber voltage (у) 
to the proper value in order to achieve a linear response due to 
neutron flux alone. This leaves room for error in the linearity 
of i.(t) if uc is not properly adjusted. Fortunately, as the 


neutron flux level increases the current contribution to i (t) 





12 - , 
Training Manuel, Instrumentation, p. 5. 
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due to gamma flux alone becomes negligible and i (t) approaches 
a linear dependance on neutron flux alone. 

Experimental results reveal that improper adjustment of 
the compensating voltage at low reactor flux levels (corresponding 
to reactor output power levels below ten watts) causes the output 
current of the ion chamber (i (t)) to become nonlinear. Figure 
2.4 shows the maximum nonlinear effects at low power levels when 
the compensating voltage is adjusted to zero volts. 

Since a linear relationship between ion chamber current 
(i (t)) and neutron flux (also reactor output power level) is 
necessary, it is best to rely only upon output currents which 
correspond to reactor power levels greater than ten watts. 

Specifications of the ion chamber are given in Appendix E. 
Since the output resistance is greater than 1013 ohms, the chamber 


may be considered as a current source. 
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Figure 2.4 Graph of ion chamber output current as a function of 


output power level for the Michigan State University 


TRIGA Mark I nuclear reactor 






СНАРТЕК ІТІ 


DESIGN OF THE REACTIVITY MEASURING DEVICE 


3.1 Formulation of the Problem 


Having presented the equations defining reactivity and its 
relationship to the output current of the ion chamber it now re- 
mains to develop a means by which р сап be measured. 

Equation (1.1) expresses the relationship between Ө 


and I. where 


С/е, 
The objective is to process this current so as to produce a voltage 
proportional to the reactor stable period. If this can be accomplished 
effectively then 95 may be related to p as in Equation (1.2) 


since 


9 
Бей 5 (1.2) 


Because this relationship is slightly nonlinear as the 
graph in Figure 2.2 illustrates, the value of p corresponding 


to a given өр пау be taken from this graph and used to calibrate 


the reactivity meter. 


15 





16 


l 
3.2 Possible Methods of Solution j 


Given the desired output р and the input i (t) several 
methods for determining р will now be considered. 

One method makes use of the circuit in Figure 3.1 where a 
d-c voltage proportional to LA is obtained from i (t) us ing 


a logarithmic amplifier followed by a differentiator and filter. 


i (t) 
5 


€ 


v4 (t) 


Amplifier tiator Filter 


Differen- ма Low-Pass 


Logarithmic 





llo 


Figure 3.1 Block diagram of a circuit which provides a d-c 
voltage proportional to the inverse of the 
reactor stable period 


In a circuit of this type the logarithm of the ion chamber 
current is first taken such that 


t/0 


At 
MS i — = р = > ---- 
v1 C) А log,gi.(t) А log, gIge А logs oly а log, ge (З) 


The voltage v4 (t) is a linear ramp function whose slope 
is inversely proportional to the negative inverse of NS Dif- 


ferentiation of this ramp voltage yields a d-c voltage V, where 


-Kd -Kd At KA 
= IM = NR == = 
V t v(t) a ( log, ofo à log, ge) lo 


3.2 
2 а J 510° 052) 
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The first two methods presented were suggested by Dr. L. Giacoletto 


of Michigan State University. The third method was developed by 
the author. 
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The result is a d-c voltage inversely proportional to 9 .. This 
voltage is then passed through a low-pass filter to eliminate any 
low-frequency noise components passed by the differentiator. 
Finally, due to the slight nonlinear relationship between 
o, and p the voltage Vm (Figure 3.1) is calibrated on a 
meter via Figure 3.6 to indicate exact values of reactivity. 
A second method for obtaining a d-c voltage proportional 


to o, involves the use of an integrator and a quarter-square 


multiplier. The circuit is pictured in Figure 3.2. 








Quarter- 
Square 


Multiplier 


Figure 3.2 Circuit for obtaining a d-c voltage proportional 


to Ө 
Р 


In this circuit the input current (i (t)) is first con- 


verted to a voltage 


t/Q 


CLE M Е р (3.3) 


using a buffer amplifier whose input resistance is much greater 


than R. 
v 


This voltage is then integrated to obtain v, (t) where 


СЕ “IR Ө E78 
VN e вс 176946 = CON Ge (3.4) 
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The quarter-square multiplier is next used to divide out 


the exponential term in v, (t) such that 





SES -IRỌ in -9 
out у. (Е) ЕС t/8 RC 
1 р 
Т Ке 
Ov 


This resultant d-c voltage is proportional to 6 and 
can therefore be calibrated on a meter using the graph of Figure 
3.6 to represent the corresponding values of reactivity. 

A third method for finding 9, involves the integration 
of a voltage proportional to the change of the input current. 


Figure 3.3 is a simplified diagram of the circuit. 


v (t) 
c 





Figure 3.3 Circuit for obtaining a discrete voltage propor- 
tional to the slope of i (t) 


The input current is first converted to a voltage MESS 


such that 


t/g 
v4 CD = i (t)R, = Т (В. е р (3.3) 


This voltage is then fed into a window comparator whose "GO" out- 
put (Figure 3.3) acts as a current sink for Y, < v4 (t) < v and 


as a large resistance whenever v(t) ZN or v4 CO > um The 


4, 


"GO" condition will allow the current i (t) to flow in the coil 
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of the reed switch and close its contacts. When this occurs the 


d-c voltage V is switched to the input of the integrator and 


T 
is integrated for the period of time that v(t) remains within 
the bounds of E and Це This period of time is proportional 


to 6 and is also proportional to the final integrator output 


voltage. The relationship between уа И 


and 8 can be found 
u P 


by starting with the output voltage integral 
Ө 


EI 
V_dt (3.6) 


Vout © RC it 


The total integration time (8,) is that period of time that 


V, < MATS < p and may be calculated by substituting m and 


V, into the equation for v(t) where 
9,/0 
ou t Фр 
У1(9,) um Rite (3.7) 
апа 
9/9, 
у109,) IR oc CE) 


Ө, апа 9, are the times for which v4 (t) is equal to and 
M. respectively. 

Equations (3.7) and (3.8) are then solved for the time 
increment 9, where Ө, = 0 - 0 апа 


u 4, 


V V 
U— an 6—) (3.9) 


lo ^ NI 


Using Equation (3.9) for the value of 2. іп the integral 





= 8 = = 
е 9 9 (4n 


of Equation (3.6), the voltage VE may be written as 


V V 
-1 u f 3 
- --- --<- у= zr 
ы eR кы 
у 0 vO 
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which is a d-c voltage proportional to the stable reactor period. 
As in the first two methods this voltage may be calibrated 
on a meter using the nonlinear graphical relationship of Figure 3.6 
to indicate the corresponding value of reactivity. 
The above method for measurement of reactivity differs 
from the preceeding two in that the output voltage is not a con- 
tinuous representation of 9, but rather a discrete value obtained 
after integrating a d-c voltage over a discrete time period 8, 
After consideration of the three alternatives the method 
chosen to measure p was the first one described using a 
logarithmic amplifier followed by a differentiator and filter 


(Figure 3.1). 


3.3 The Circuit and Theory of Operation 


With the reactivity measurement solution method chosen, 
the circuit of Figure 3.4 was designed to produce a d-c voltage 
proportional to 95 

The transfer function of the circuit is 


52615) 1 


V = Y (o татын o cL 
out 1 (к1С‚в+1) (К/С„5+1) (В С 531) 


(3.11) 


where V4 (s) is the frequency domain expression for vı (t) given 


as 


' -А log 9 (15/1, ) А Тов. (е 

У (5) =т= Наше этсе = 
S ane 
P 


(3.12) 


and A is the logarithmic amplifier voltage constant determined 


approximately by 
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К + 5000 + R -3 
26 X ) 
А = — MÀ R Al 2 S Volts (3.13) 
3 . 


where Ra is a temperature dependant resistor within the logarithmic 


amplifier ca 

I, is the log amplifier reference current determined 
primarily by RUE 

The concept by which the circuit works is as follows: 
Using the ion chamber current as the input, the logarithm of i (t) 


is first taken. The output voltage (vi(t)) is a linear ramp 


function since 


iagt) I t/g 
v (t) = -A log, 9 PESE -А log, ¿e P) 
k k 
10 At 
= -À 108, Те” e, 28104 (3.14) 


v] Ct) is then differentiated to produce a d-c voltage V, л 


| Modular Logarithmic Amplifier Model 4116 (Burr-Brown Research 
ee”, 1971), р.3. 


E A more complete discription of the logarithmic amplifier and 


its parameters is given in Appendix A. 


16 Р . , i , 
Since the input current (i, (t)) is an exponentially increasing 


function and since the circuit (Figure 3.4) used to process this 
function has a limit to the absolute magnitude of the input, i (t) 
can only remain "on" for a discrete period of time E Therefore, 
in order to be absolutely correct when defining i the gate 


function notation should be used where 


o = 1 Ostsee, 


G Е - = 
= 9/2) 0 с<0, с> 0. 


Ө, is the time of the input pulse. 
Using this notation i (Œ) becomes 
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I 
d d 0 At 
A O o ЖИ: 18106) 
R.C,A 
241 
— 3215 
a, 108, 9° Ga 


This d-c voltage is then filtered to eliminate any low- 
frequency noise enhanced by the differentiator. The resultant 
voltage уа is inversely proportional to the reactor stable 


period and approximately proportional to the reactivity such that 


КСА M КСА 





a = ЕСЕ T log, ge (3216) 

As previously discussed (Section 2.2) the relationship 
between 9 and p is nonlinear. Calibration of the meter may 
be accomplished by producing known reactor periods and marking 
the meter deflection with the corresponding values of p via the 
graph in Figure 3.6. 


A comprehensive analysis of the circuit design is given in 


Appendices A, B, C and D. 


3.4 Calibration 


The output voltage in the circuit of Figure 3.4 is a d-c 


voltage proportional to the inverse of en and approximately 


t/s 
EE E, (50.17) 


If the gate function notation is carried with i (t) two finite 
strength impulses will appear when v(t) is differentiated in Equa- 


tion (3.15) due to differentiation of the gate function. 
Although these finite strength impulses are of importance (partic- 
ularly since they show up on the output when i (t) is switched on 


and off) they shall not be carried in the analysis. ТЕ should be under- 
stood, however, that i (t) is on for a finite period of time. This 


period of time (0,) is discussed further in Appendix B. 
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proportional to p where 


КЕС 
€ BUE log ESTO ү. 
10 8 


3 
М ТЕ 9, 108, OF ( 2 


As previously mentioned the relationship between ЫН апа 
p is nonlinear. Figure 3.6 illustrates this relationship for the 
Michigan State University TRIGA Mark I reactor. 

The range of reactor stable periods for which the meter is 
to be calibrated is 10 sec < a s 40 sec. This range corresponds 
to a range of reactivity between 40 and 19 cents respectively. 
(Figure 3.6) 

Calibration of the reactivity meter is accomplished with 


the apparatus in Figure 3.5. 






Reactor Power 
Level Indicator PR Reactivity Meter 


Ion Chamber 
Figure 3.5 Apparatus used to calibrate the reactivity meter 


The reactor is first made critical with the ion chamber 


current equal to Lo: A control rod is then pulled out a discrete 


distance and the ion chamber current becomes 


С/е, 
i (t) = 12 (ШШ) 





25 


Тһе corresponding d-c voltage on the reactivity meter will 


be 





log (3.16) 


10° 
This voltage is then marked on the meter. 

In order to mark the reactivity corresponding to Noe 
the graph of Figure 3.6 is used. As this graph indicates the 
reactor stable period must be known. This value is determined by 
using the reactor power level indicator and a timer to find the 


power level doubling time JE The period can be calculated from 


the doubling time since 





ox (3.17) 


This process is repeated until the range of p is 


adequately covered. 





p, Reactivity (cents) 
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О? Period (seconds) 


Figure 3.6 Graph showing the relationship between reactivity and 
reactor stable period for the Michigan State University 
TRIGA Mark I nuclear reactor (Courtesy of Gulf General 
Atomic) 





СНАРТЕК ТУ 


CONSTRUCTION OF THE REACTIVITY METER 


4.1 Circuit Tests, Performance and Specifications 


The circuit of Figure 3.4 was constructed and tested. 
Figure 4.1 shows several photographs of the instrument. 

Testing of the circuit was accomplished by connect ing 
the meter input to the ion chamber via a coaxial cable. The re- 
actor was then made critical and a control rod withdrawn to pro- 
duce a reactor period within the desired range. The range of 
the meter was checked for consistancy and variation due to noise. 
These tests revealed the meter deflection to be approximately 77 
percent for a 10 second period and approximately 24 percent for a 
40 second period thus adequately covering the desired range of 
Qe: The maximum variation in meter reading due to noise was less 
than + 0.4 percent of full scale. 

The minimum reactor power level for which the circuit was 
designed to operate is ten watts. This corresponds to an ion 
chamber current, Ig < ЗаА. Values of input current below this 
level will put the logarithmic amplifier in saturation and there- 


fore cause Ye to register inaccurately on the meter. 


ut 
The upper power level was chosen to be 200 watts resulting 


in an ion chamber current of 60 nA. This value is not an absolute 


upper limit and the meter will continue to read accurately up to 
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a reactor power level of 1000 watts (i (t) = 280 na) at which time 
the negative temperature coefficient inherent in the reactor will 
cause the meter to lose accuracy. 

The final specifications of the reactivity meter are listed 


in Table 4.1. 


Table 4.1 Specifications of the reactivity meter 


E 
e 9% 
0 


10 seconds 


Input: Exponentially increasing current «+. i (t) zx 


8 ә ә ө ее ә  ө е ә е ә ә е ә ә е ә ә ә ә ә ә ә ә е ә ә ә ә ә ә ә е е өе ө ө э ә 40 seconds 
p max 


a iiie. v RE S Ка А БАСА DOR Amperes 


5 тіп 


60 X 1054 Amperes 


.44...ш.шш.вшш................- ГОЯЗСЕТУТЕУ (р) in cents 
еее ә ә ә ә е ә е ә ә ә ә ө ә ә ә ә е ә Өө 900.00.» 19 cents 


40 cents 


. * 0.5 percent 





Since the values of р оп the meter are calibrated from 
a graph, the only error will be that variation caused by noise on 


the output voltage of the filter (Figure 3.4). 






CHAPTER V 


CONCLUSIONS 


The control rod calibration curve in Figure l.l is of the 
type for which the reactivity meter was designed to produced. 
Original design specifications called for the measurement of periods 
between 10 and 40 seconds corresponding to reactivities between 
40 and 19 cents. As the graph indicates, these values cover only 
a small portion of the total rod calibration curve. 

In order to cover the entire range of the calibration 
curve the measured values of rod reactivity must be added in incre- 
ments within the range of the meter. For example: Reference to 
.the graph of Figure 1.1 shows the reactivity corresponding to a 
rod position of 300 to be approximately 60 cents. If at this 
position the reactor were made critical and the rod pulled out to 
a position of 350, a reactivity of about 40 cents would be recorded 
on the meter. In order to get the total rod worth at this new 
position this value of reactivity would have to be added to the 
previous total reactivity of 60 cents and would therefore yield 
the new total reactivity of approximately 1 dollar. 

It should also be noted that since the total rod worth is 
determined by the sum of incremental reactivities the error in 
measurement of each incremental rod worth becomes very important 


and the total error increases with each sum. 
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The overall results of the tests made on the instrument 
reveal that the reactivity meter works well for the range of input 


currents and reactor periods for which is was designed. 
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APPENDIX A 


DESIGN OF THE LOGARITHMIC AMPLIFIER! 


The logarithmic amplifier used in constructing the circuit 
of Figure 3.4 is a Burr-Brown model 4116 modular logarithmic 
amplifier. Its purpose in the circuit is to take the logarithm 
of the input current i (t) where 

</9 


i Ors ie P (1.1) 


The resultant output is a linear ramp function given by 


i (t) Т, СЯ 
v(t) = -A 108, TT i -А 10840 —- e 
k k 
_ 10 Ar 
= -A 108, 9 E - 8, log, ge (3.14) 


There are two design parameters for this amplifier. These 


are the voltage scale factor A and the reference current I. 


Figure A.1 shows the logarithmic amplifier with exterior connections. 


The electrical specifications are given in Table A.1. 


E The design equations, graphs, figures and specifications in this 


section were taken from the pamphlet Modular Logarithmic Amplifier 
Model 4116 (Burr-Brown Research Corp., 1971), pp. 1-5. 


33 









34 


Table A.1 Electrical specifications of the Burr-Brown modular 
logarithmic amplifier model 4116 (current input only) 


Accuracy 
Accuracy, percent of full scale with current 

source input for 0.4nA < І. < 400 PAE T + 1 percent 
Input 
Current SOUmeemmen DUC e655 665 a old uence one 400 pA to 400 yA 
Absolute Maximum Current ......... ос ОЕ cess & X0 mA 
Reference Current Range (Ij). eee sa . 400 nAs I s 40 „A 


Voltage Scale Factor Range (A) ............ 2/3 V < А < 10 V 


Output 
биттеп, ee a к a. КЕ NE... е аси + mA 
a РТР + 10 У 
Output Impedance at A =D man. ee ee 10 Ohms 
Stability 
Scale Factor Drift (AA/°C) ...... cos ....... +`0.0005 А/С 
Reference Current Drift (A1, /"C) 

poU EE Wer caes + 0.003 1, /%с 


Input Offset Current Drift о pA at 2556. l рА/% 


Input Noise - Current Input .... 1 pA rms, 10 Hz. to 10k Hz. 


Power Supply 
Rated Voltage sen... snes ee en 5, Е СОСЕ 


Supply Drain (Quiescent) e a a 20550020 01 + 22 шА 
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R 

X 
++ 15 V 

3 ts R 

Current Input да Х ВВ - 4116 

------ 5 
Voltage Input D Log Amplifier 
, 4 D V 
(optional) ә 1 


+15V6 dec offset 5b- 15 V 
adjust 


Figure A.l Burr-Brown modular logarithmic amplifier model 
4116 with external connections 


A design that takes advantage of the full range of out- 
put voltage is best. The values of Ro and ES required to 
accomplish this may be obtained by first solving Equation (3.14) 


tfemmthe full + 10 volt range where 








T | 
V4 ге - -А lo8, эш = 10 Volts (A.1) 
k 
T 
Лл да дағ x = = -10 Volts (А.2) 
k 


These equations may be added to solve for I. with the 


result 


^^ 


Е ZU 
к, SE. Pin ( ) 


The corresponding voltage scale factor (A) may next be 
found by substitution of Equation (A.3) into Equation (A.1) for 


Г This provides an expression for А where 






36 


ay ашык (A.4) 


І 
logig max 


Los 
min 





In Chapter IV it was established that the value of Lum 
would be 3пА апа I would be 60nA. Substitution of these 


values into the expressions for I and A yields 


I. = 1229648 (A.5) 


A = 15.75 Volts (A.6) 


At this point a problem is encountered. Figure A.2 shows 
the graphical relationships between I, А and their design re- 
sistors. Reference to these graphs and to the specifications in 
Table A.l indicate that the values of I. and A in Equations 
(А.5) and (A.6) are outside the design limits. 

A compromise must be made to reduce the output voltage 
range in order to put I. and A within their respective tolerances. 
For the final design a value of R_ = 11.3 MQ was selected. This 


corresponds to a reference current of .56,A. 


Using this value for I the best value for A can be 





k 
calculated by letting the output be 10 volts for - E Sub - 
stitution into Equation (A.l) yields 
A M LC MN NE (A.7) 
Ro. -9 
log min ee 3x10 
PME MEINE 


Reference to the graph in Figure A.2b indicates a value of 
28 Ki for Ба will make А equal to 4.5 volts. 
The actual output range of the logarithmic amplifier may 


now be calculated where 






100 yA 
10 pA 
100 nA 
10 nA 
10 k 100 k IM 10 M 100 M 
R (ohms) 
У 
(а) 
10 V 
8 V 
6 V 
< 
4 V 
Zo 
0 
0 20 k 40 k 60 k 80 k 
R (ohms) 
X 
(b) 


Figure A.2 Design curves for the Burr-Brown model 4116 
log amplifier, (a) Reference current vs. R , 
(b) Voltage constant vs. R. У 
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ee = 10 Volts (А.8) 
it -9 
60x10 
. = -A log EX = _4.5 log. . z——— с 
] min 10 In 10 56x10 6 
= 4.36 Volts (A. 9) 


The design parameters of the logarithmic amplifier may 


be summarized where 


A = 4.5 Volts 
R = 28. k Ohms 
y 
Ik = лб 
R. = 11.3 M Ohms 
Input Range (i (t)) 3nA to 60nA 


Output Range (v4 CO) 10V to 4.36V 
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APPENDIX B 


DETERMINATION OF THE RAMP FREQUENCY COMPONENTS 


The ideal output of the logarithmic amplifier in the circuit 
of Figure 3.4 is a linear ramp function of the form 
I 
Zo _ At 


А 
I J log} o£ (3.14) 


v} (t) = -À logi 


ж 


Differentiation of this function will yield a d-c voltage 





where 
d d E At 
у, - “RC, DICE) = “RC О log. o I LE log, ge) 
P 
КСА 
БЕ: 2 [| 
= log, ge (3 23) 


9 
In order to obtain the true differential as in Equation 
(3.15) it is necessary to differentiate all frequencies of the 
ramp function v(t). This therefore requires the knowledge of 
the range of ramp frequency components. 
The highest frequency components occur when the ramp func- 
tion is on for the shortest period of time. This corresponds to 
a reactor stable period of ten seconds. 
The "оп" time for this ramp may be found by substituting 


the limits of i (t) into Equation (1.1) and solving for 9, where 


D. pM PUE (1.1) 


E 





&0 


от 
o, /10 
ба ee: 
Ө, = 10/,п 20 = 30 seconds 
The Fourier series half-range cosine expansion for 
f(x) = x as a recurrent triangular wave is 
k 4k TTX 1 TX 1 5 
= = == — + -- —— + -- Бадаа 
f(x) x > > (cos Е g COS = 25 COS = 
TT 
2n-1) E 
( Thc N (В.1) 


Substituting the ramp function of Equation (3.14) for 


f(x) the Fourier series expansion for vi CO is 


2 І Ө А 1ор, е 
Ас qr c E 10 
(С) = «<A log ше log, e = -A log mA uu uc жаланы 
2 
1 10 Іі, 9, 10 10 I, 95 m 
49 A 
: log, e|cos ЗЕ + er Se 
- 10 Өнө 9 GNE 25 я 
па. t 4 
p min 
1 7gt , 1 Ont , 1 Lint 
Stale COG Uit CO съл ес tC OG 
49 ОЕР u, il o. 
coc = 2 cos T me Оса Ена O, = 30 sec. 
(2n-1) t E S ae 


The radian frequency components for the minimum duration 
ramp and the respective amplitudes of each (Р) can be determined 


from Equation (B.2). The values are 


1 
$ Samuel M. Selby, CRC Standard Mathematical Tables (Cleveland, 
Ohio: The Chemical Rubber Co., 1960), p. 147. 
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49,4 log, ge 


в 


Q4 = = 0.1048 гаа/зес P. = ——————— 27.9 
1 » 5 
9t | т Ө . 
p min 
49 A log, ,e 
шз = a = 0.314 rad/sec, Py = E = 0.878 
t 9п 9 i 
p min 
49 A log, е 
P E RA O оү 
E Q 5 2 
t 25т 8 . 
p min 
49 A loge 
шу = шаа 2 = — AL = 0.161 
9c 491 0 ) 
р шіп 
48 À 10р. е 
Wo и - 0.943 rad/sec, Py - = 0.0975 
t 81п 9 ; 
p min 
48 A log, „е 
шуу = Inn Mud SR Poe A = 0.0653 
9. win o. 
p min 


Since the amplitudes of higher frequency components are 


less than 1/144 of the amplitude of ©, it is reasonable to omit 


them from the analysis. The value of w is therefore the highest 


11 


radian frequency to be differentiated. 
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APPENDIX С 


DESIGN OF THE DIFFERENTIATOR 


Figure C.1 is a magnitude plot for the real and imaginary 


components of an ideal differentiator whose transfer function is 


uu) 
= Zu O (C.L) 
M in (w) 
Se а (w) y vous 
ч ur (yu a 
slope = 
(a) (b) 


Figure C.1 Transfer function plot for an ideal differentiator, 
(a) Imaginary part, (b) Real part 


A differentiator circuit with the type of transfer function 
in Equation (C.1) has the disadvantage of enhancing noise. 

If the range of frequencies to be differentiated is limited, 
it is possible to reduce unwanted noise by using a compensated 


differentiator whose transfer function is 


1 
2 Albert S. Jackson, Analog Computation (New York: McGraw-Hill 
ВооК Со., 1960), р. 147. 
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V 
out іш jKw (C .2) 
vo GKo + 1) GKo + 1) | 
where 
V (w) 2, 2 
2K 
Re | -2 — | = — U (C .3) 
Vin% ie 
and 
ml. eu? | _ ma - Ku) T 
= 5 ; 
Vin @ KDE T AK S. 
These real and imaginary components are plotted in Figure 
6222 
V 
Т E ра RER 
C A us 
| | 
u) 


| (a) (b) 


Figure C.2 Transfer function plot for a compensated 
differentiator, (a) Imaginary part, (b) Real part 


In Appendix B it was determined that the maximum radian 
frequency to be differentiated in the linear ramp function was 
шуу = 1.15 rad/sec. It is therefore desirable to make 1/K greater 
than 1.15 rad/sec. 

Figure C.3 is a circuit diagram of the differentiator. 


20 y 
КСА Linear Integrated Circuits (RCA, 1967), pp. 75-78. 
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Figure C.3 Equivalent circuit of a compensated differentiator 


The transfer function may be derived as follows: 





Define: 
Ry 
Ze 7 R63 TORO, * 1 $ 
22 
1 
Z =R, + 
js 1 ЗС € 
1 
Уз» m and Vol may be written as 
ZORIZ и Т 
ү = ү f oi” 1 “a O (С 
TEE T + Ши + + 7 
е іп 2. (2. 2094) і Ze zi NIA 
Vi 00 (С 
= Г V - C 
Vg, = A (0) (V, - VD 
Substitution of Мы, апа Уш into the expression for 
Vol yields 
У. А р 
V Б E o m ET (С 
4 + + + + A 
01 (2 ғ 21221 ZZ, Loi Z) 0 w)Z 2, 


5) 


.6) 


.7) 


.8) 


2) 


.10) 


If R, > œ the exact expression for the transfer function 


4, 


сап be written as 
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A 7 
Bus ZA o «927, 


" + + Y 
Vin (шуға: шасы астра 


(С.11) 


This expression may be simplified if the following assump- 
tions are made: 


Assume: 


2, >> Z Sz (C212) 
Zoi << 2, (15) 
Ze << Âg (w) (С.14) 
2. << А.(ш) (C .15) 


2 
Equation (C.9) may then be written | 


(С.16) 


When the expressions for 2. апа Zr are substituted into 


this equation the differentiator transfer function becomes 


x : -jJR4C W 
үшы аты ЫЗЫ с 
їп Ju 17295 


If Ric, = RC, then the magnitude plot in Figure C.2 


will correspond to Equation (C.17) with 1/К = 1/R,C,- 
For a differentiator error of one percent maximum it is 


necessary to establish a corner frequency about ten times the 


maximum frequency to be differentiated. In other words 


1 1 
——— ш ----- - 10 wu = 10 uv (C.18) 
К.С. RC, max 11 


21 
Ibid., pp. 75-78. 
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22 


where = 1.15 rad/sec and was determined in Appendix B. 


zn 
A trade off of accuracy vs. noise is encountered if Equa- 
tion (C.18) is used to design the differentiator. The prupose of 

the low-pass filter in Figure 3.4 is to attenuate the noise passed 

by the differentiator. It was discovered that for 1/R4C, = 10 41 

= 11.5 rad/sec the filter could not effectively reduce the noise 

without causing an excessive delay in the d-c response due to its 

long time constant. A decision was made to limit the filter time 
constant to 0.6 second. This meant the corner frequency of the 

differentiator had to be reduced until the noise level was tolerable. 

The experimental resultant corner frequency was 1/R,C, = 3.8 rad/sec 
or about three times as great as the highest frequency to be dif- 
ferentiated. 

The magnitude of the error caused by inexact differentiation 


of the maximum ramp frequency will be small since the amplitude of 


this frequency component is only 1/121 of the major ramp frequency. 


22 
Jackson, Analog Computation, p. 147. 
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APPENDIX D 


DESIGN OF THE FILTER 


Figure D.1 is a schematic of the output filter in Figure 
с 
Тһе filter is a simple RC low-pass filter with a unity 


gain output buffer stage to prevent loading. 





d-c offset 
| adjust 
= - 15 V 


Figure D.1 Output filter for the reactivity meter 


The transfer function for the filter is 


V 
Td = 1 (Ф.1) 
Vig (D) JR £C 20 +l 
Since the output voltage of the differentiator (V,) is 


a d-c voltage with superimposed low-frequency noise it is desirable 


to make the filter time constant long enough to reduce this noise 
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while at the same time have a reasonable response time to the d-c 
voltage. 

Final tests of the reactivity meter indicated that a filter 
time constant of RC. = 0.6 seconds reduced the output noise level 
on the meter to less than + 0.4 percent full scale. This value 


ГОТ RC. also gave a d-c response time of about 3 seconds. 
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APPENDIX E 


SPECIFICATIONS OF THE COMPENSATED ION CHAMBER 


The TRIGA Mark I nuclear reactor at Michigan State University 
contains a Westinghouse type WL-8501 compensated ion chamber. Table 
23 
E.1 tabulates the specifications of this chamber. 


Table E.1 Specifications of the Westinghouse type WL-8501 
Compensated ion chamber 


Maximum Voltage Between Electrodes (dc) 

Operating Voltage for Boron Chamber (dc) UU 

Operating Voltage for Gamma Chamber (dc)........ 0 бо -37 V 
l3 


Output Resistance, Signal to Case, Minimum ........1O Ohms 


Output Capacitance, Signal to Case, Approximately 


Thermal Neutron Flux Range 2 5x10? to 1.5x10l enc / sec 
1 


Thermal Neutron Sensitivity КИ тото ; А/ пу 


Gamma Sensitivity, Conan or A/R/hr 





Since the output resistance of the ion chamber is greater 


13 
than 10 ° Ohms the chamber may be treated as a current source. 


Westinghouse WL-8501 (New York: Westinghouse Electric Co., 
1965), pp. 1-2. 
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